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ABSTRACT 


Integral  equations  are  derived  to  model  the  solid 
surfaces  and  the  inlet  condition  for  both  external 
and  internal  flows  with  an  axially  rotating  element. 
Vortex  filaments  are  used  to  simulate  the  lifting 
surfaces.  Global  iterations  for  surface-panel 
source  densities  and  vortex  strengths  are  used  in 
conjunction  with  the  Neumann  iterations  for  solving 
integral  equations.  A  special  kernel  treatment  is 
adopted  in  the  solution  scheme  of  the  Integral 
equation.  Solutions  converge  within  10-20  iterations. 
Calculated  results  and  measured  data  for  a  propeller 
tested  in  open  water  are  in  good  agreement.  Predic¬ 
tions  are  also  made  for  an  axial-flow  pump.  The 
overall  performance  parameters  for  the  pump  agree 
within  a  few  percent  of  the  design  values. 


ADMINISTRATIVE  INFORMATION 

This  research  was  funded  by  the  Surface  Ship  and  Submarine  Technology  Program, 
PE62121N  Block  NDIA  and  PE62323N  Block  ND3A,  respectively,  administered  by  the  David 
Taylor  Research  Center,  Code  0124. 


INTRODUCTION 

It  has  become  a  requirement  in  the  design  process  to  obtain  numerical  solutions 
for  an  external  or  an  internal  flow  past  a  rotating  element  with  lift-generating 
surfaces;  e.g. ,  marine  or  aircraft  propellers,  turbine  or  pump  Impellers.  Although 
a  complete  viscous  solution  technique  would  be  a  most  desirable  tool  to  achieve 
design  accuracy,  a  number  of  difficulties  exist  in  its  applications.  These 
difficulties  Include  computer  time  and  memory  requirements,  numerical  solution 
procedures,  turbulence  modeling  and  lack  of  experience  in  handling  such  complex 
computer  codes.  Because  of  its  modest  complexity  and  reasonable  computational  time, 
an  inviscid  solution  is  therefore  often  pursued  in  the  direct  design  process.  In 
many  cases,  the  accuracy  of  the  Inviscid  solution  is  adequate  for  design.  For  the 
inviscid-f low  model,  Miranda  [l]  has  compared  the  pros  and  cons  between  the  potential 
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and  the  Euler  solution  procedures.  In  general,  the  selection  of  a  numerical  solution 
procedure  Is  often  case  dependent,  and  no  conclusive  argument  can  easily  be  made. 
Significant  reductions  of  the  efforts  In  grid  generation,  engineering  and  computer 
time,  and  computer  memory  requirements  promote  the  adoption  of  the  potential 
solution  using  a  surface-panel  singularity  method  [2]  instead  of  solutions  using 
either  a  finite  difference  or  a  finite-element  procedure  [3]. 

The  vortex-lattice  method  [4]  has  been  used  to  model  the  axial  flow  through  a 
marine  propeller  for  many  years.  Realizing  the  shortcoming  of  this  approach  In  rep¬ 
resenting  arbitrary,  thick  surfaces,  reference  [4]  was  extended  [5,6]  to  Include  the 
Influence  of  the  dominant  thickness-effect  components;  e.g.  hub  and  duct  for  the 
propeller.  In  recent  years,  Hess  [7]  has  successfully  applied  his  earlier  develop¬ 
ment  [8]  to  aircraft  and  marine  propellers,  using  a  source  density  on  all  on-body 
panels  and  a  dipole  density  on  lifting-surface  panels.  Maskew  [9]  extended 
Morino's  [10]  approach,  which  uses  Green's  Identity  to  determine  a  unique  source- 
doublet  specification,  to  the  helicopter  rotor  problem.  Nathman  [ 1 1 ]  also  applies 
Morino's  approach  to  the  internal-flow  problem.  The  drawback  for  Nathman' s  applica¬ 
tion  is  that  Green's  Identity  is  only  valid  for  a  closed  domain.  This  restriction 
imposes  a  requirement  of  boundary  conditions  at  the  downstream  outlet  section  for 
the  internal-flow  calculation.  In  general,  the  outlet  conditions  are  part  of  the 
solution,  and  it  is  difficult  to  predetermine  them. 

In  this  report,  integral  equations  for  surface-source  densities  are  derived  and 
compared  for  both  external  and  internal  flows  with  an  axially  rotating  element.  For 
internal  flows,  integral  equations  are  formulated  for  nonuniforra  inlet  flow  and 
"free"  (i.e.,  no  prerequisite  conditions)  outlet  flow.  Various  strength  distribu¬ 
tions  for  the  vortex  filaments  are  examined  and  compared  to  achieve  an  optimal  rep¬ 
resentation  of  the  lifting  surfaces.  A  numerical  formula  using  the  Blot-Savart  law 


for  the  evaluation  of  the  velocity  field  due  to  a  vortex  segment  Is  derived.  Global 
Iteration  between  the  source  and  the  vortex  strengths  Is  used  in  conjunction  with 
the  Iterative  scheme  for  solving  the  Integral  equation.  A  special  treatment  [ 1 2 ]  to 
eliminate  the  singularity  of  the  kernel  of  the  integral  equation  Is  adopted.  Two 
numerical  examples,  one  for  a  marine  propeller  and  one  for  an  axlal-flow  pump,  are 
presented. 


MATHEMATICAL  MODEL 

STATEMENT  OF  PROBLEM 

The  problem  being  considered  is  the  determination  of  the  incompressible  axial 
flow  field  through  a  rotating  element  with  lifting  surfaces.  The  flow  can  be  either 
unbounded  or  bounded  by  a  duct.  The  domain  of  Interest  V,  see  Fig.  1,  may  or  may 


BOUNDED  S  *  Sr  +  Sc  +  Sj 
Fig.  1.  Definition  sketch. 

not  be  simply  connected.  The  boundary  surfaces,  S,  may  consist  of  rotating  or 
nonrotating  inner  surfaces  Sr  and/or  outer  boundary  surfaces  and  Sc:  S  =  Sr  +  Sc 

+  Sj.  The  boundary  conditions  at  the  outer  surface  Sc  and  the  nonrotating  inner 
surfaces  Sr  are  flow  impermeable  conditions,  l.e. , 
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where  N  is  in  the  surface  normal  direction  and  directed  into  V-.  For  the  rotating 
inner  surface,  the  fluid  velocity  in  the  N-direction  is 

VN  =  (to  x  r)N  (lb) 

where  <o  is  the  rotational  velocity  and  r  is  the  positional  vector.  The  outer  bound¬ 
ary  surface  Sp  is  a  permeable  surface  used  to  describe  the  nonunif orraity  of  the 
inlet  condition  for  internal  flow, 

VN  =  vin  .  (lc) 


METHOD  OF  APPROACH 

The  method  of  singularities  is  utilized  for  the  solution  of  this  rotat ing-f low 
field  problem.  To  satisfy  the  boundary  conditions  (1),  sources/sinks  are  distrib¬ 
uted  on  all  the  inner  surfaces  and  also  on  the  outer  surfaces  for  a  bounded  flow 
case.  To  satisfy  the  wake  tangency  condition  at  the  sharp  trailing  edge  of  the 
lifting-surface  element,  as  shown  in  Fig.  2,  vortex  filaments  are  distributed  along 
the  mean  lines  of  the  cross  sections  of  the  lifting  surfaces  as  well  as  in  the  wake. 


REPRESENTATIVE 
U-SHAPE  VORTEX  FILAMENT  r„ 

V  TRAILING 

^MEANUNE  yORTE* 


SOUND 

VORTEX 


7  &  :  '  4-V— 


Fig.  2.  Vortex  distribution  on  a  lifting  surface. 


Denote  the  velocity  vector  due  to  the  surface-distributed  sources  by  Vfj  and  that  due 


to  all  the  vortex  filaments  by  Vp.  We  then  have 
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(2) 


V  -  V<J>  =  vM  +  vr  , 


and 


vM  -  v/s  -  ^  dSq  , 
rpq 


rij(Q)  x  rPQ 

vr  -  ?  U - • 


i  1  L 


4ir 


PQ 


(3) 


(4) 


where  <j>  is  a  scalar  potential  for  the  resultant  potential-flow  field,  M(Q)  the 
strength  of  the  source,  the  distance  between  a  field  point  P  and  a  source/vortex 

point  Q  on  S  or  on  mean  lines  of  the  lifting  surface  cross  section,  and  r i j (Q )  the 
strength  of  each  U-shaped  vortex  filament  at  the  ith  chordwise  strip  (total  I 
strips)  and  jth  spanwise  strip  (total  J  strips).  Both  M(Q)  and  r i j (Q )  are  unknowns. 
An  integral  equation  Is  formulated  to  determine  the  strengths  of  the  M  distribution 
such  that  boundary  condition  (1)  is  satisfied.  A  set  of  linear  equations  with  J-l 
unknowns  are  set  up  for  the  strengths  of  the  distribution  such  that  the  wake 
tangency  conditions  at  K^,  K2,...  Kj_i  (Fig.  2)  are  satisfied. 

The  flow  field  in  V-  can  therefore  be  computed  by  using  Eqs.  (2)-(4).  The  local 
pressure  coefficient  Cp>  acting  on  the  lifting  surfaces,  can  be  calculated  by  using 
the  steady-state  Bernoulli's  equation. 


1  - 


(5) 


l/2p  (r2ajZ+V002  ) 


where  V  Is  the  absolute  velocity  from  Eq.  (2),  V,*,  is  the  incoming  flow  velocity,  and 
p  Is  the  fluid  density. 

INTEGRAL  EQUATION  FOR  SOURCE-STRENGTH  DETERMINATION 

When  the  field  point  P  Is  on  S,  Eq.  (2)  can  be  written  In  the  N-directlon  as 
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(6) 


~  ac,i.(p)  +  a0(P)  +  ag(P)  (8) 

represents  the  total  flux  due  to  a  "sink"  of  unit  strength  at  P  (referring  to 
Eq.  (3))  or:  S  entering  the  region  V  through  the  remainder  of  S+S0.  When  P  is  on 
Sc+St,  then  +  a0  =  2tr  and  ug  =  0;  when  P  is  on  Sg,  then  uc>i  +  a0  =  4ir  and 

otg  =  -2ir .  When  P  is  on  Sc+S^  or  on  S^, 

/  K(Q,P)  dSQ  =  2tt  .  (9) 

S+S0 

With  the  aid  of  Eq.  (9),  Eq.  (6)  can  be  cast  in  the  form 

4tt  M(P )  =  f(P)  +  Js  [  M(Q)  K(P,Q)  -  M(P)  K(Q,P)]  dSQ  (10a) 

for  the  unbounded  flow  case  and 


-6- 


(10b) 


2*  M(P)  -  f(P)  +  Js  [M(Q)  K(P,Q)  -  M(P )  K(Q,P)]  dSQ 
+  [2*  -  a0(P)]  H(P ) 

for  the  bounded  flow  case,  where  a0(P)  is  the  solid  angle  of  S0  with  respect  to  P, 
as  defined  In  Eq.  (8),  and  f(P)  ■  Vfj(P)  -  Vpjj(P).  Landweber  [  1 2 ]  has  shown  that 
the  limit  of  the  mean  value  of  the  Integrand  In  Eq.  (10)  is  zero  when  Q  approaches  P. 
Thus  the  values  of  K(Q,Q)  need  not  be  computed. 

VORTEX  FIELD  AND  ITS  DISTRIBUTION 

The  flow  field  generated  by  a  vortex  line  can  be  evaluated  by  the  Blot-Savart 
law  as  shown  In  Eq.  (4).  For  numerical  evaluation,  the  vortex  line  Is  discretized 
Into  a  number  of  straight  line  segments  along  its  centerline,  and  each  segment  Is 
represented  by  a  vortex  filament.  Consider  one  of  the  filaments  AB  as  shown  In 
Fig.  3.  The  velocity  field  at  a  field  point  P  Induced  by  this  filament  Is  given  by 
the  Blot-Savart  law, 

B  ^  dh  x  r  (11) 

A  4x  r3 

where  r^g  Is  the  constant  strength  of  the  vortex  segment  AB,  and  dh  »  AB. 


P 


Denoting  AP  by  £,  PB  by  5  ,  and  performing  the  integration  in  Eq.  (11)  yield 


Vab  -  I 
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(12) 


VAB 

where 


rAB 

4u 


g  (hx£  ) 


g 


h2C2-(h*C)2 


( y +iei ) . 


(13) 


When  P  Is  on  the  segment  AB  or  on  its  extension,  i.e.,  either  £  =  0,  or  £ '  =  0 

2  2 _ 2  * 

or  h  ^  -(h*£)  =  0,  the  current  filament  AB  has  no  contribution,  i.e.,  V^g  =  0. 

Hence,  Eq.  (11)  is  regular  for  the  entire  domain. 


NACA  0012,  o  =  10° 


Fig.  4.  Calculated  pressure  distributions  for 
NACA  0012  foil  at  a  -  10°. 


(b) 


In  order  to  satisfy  the  wake  tangency  condition  at  the  trailing  edge  of  the 
blade  surface,  vortex  filaments  shown  in  Fig.  2  are  distributed  in  the  blade  sections 
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Cions  as  well  as  In  Che  wake.  Various  distribuCions  can  be  used  Co  achieve  this 
goal*  For  a  Chin  blade  representation,  as  commonly  used  In  Curbomachlnery  blade 
design,  a  mean-line  dlscribuCion  is  used. 

Numerical  experlmencs  were  performed  for  various  plausible  distribuCions  Co 
obtain  an  optimized  mean-line  distribution.  Two  2-D  NACA  foils,  i.e.,  NACA  0012  and 
66-209,  were  selected  for  examination.  The  NACA  0012  foil  has  zero  camber  and  was 
analyzed  at  an  angle  of  attack  a  ■  10®.  The  NACA  66209  foil  has  relatively  large 
camber  and  is  analyzed  at  a  ■  2°.  Figures  4  and  5  show  four  different  computed 
pressure  distributions  without  trailing  vortices  and  the  “exact"  inviscid  pressure 
distribution  obtained  by  the  conformal  mapping  procedure  [ 13 ] .  The  first  vortex 
distribution  is  the  simplest  selection,  i.e.,  a  single  vortex  at  one-quarter  of  the 

NACA  66209,  o  =  2° 


Fig.  5.  Calculated  pressure  distributions  for 
NACA  66-209  foil  at  a  -  2®. 
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chord  based  on  James'  [ 14]  recommendation.  Surprisingly,  the  prediction  is  reason¬ 
able,  except  near  the  vortex  location.  The  second  distribution  is  a  constant 
strength  profile  distributed  chordwise  between  10  and  85%  of  the  chord.  Because  of 
the  sudden  drop  to  zero  In  vortex  strength,  a  peak  in  Cp  Is  predicted  at  x/c  ■  0.85. 
In  order  to  eliminate  this  peak  In  Cp  associated  with  the  constant  distribution,  a 
linear  distribution,  with  unit  strength  at  x/c  **  0.10  and  zero  at  x/c  *  0.85,  was 
examined  for  the  third  distribution  (omitted  In  Fig.  5  for  clarity).  Figure  4b  shows 
not  only  a  smooth  distribution,  but  also  a  distribution  conforming  with  the  "exact" 
distribution.  The  last  distribution  in  which  the  vortex  strength  is  proportional  to 
the  local  foil  thickness,  shown  in  Figs.  4  and  5,  gives  the  best  agreement  with  the 
exact  solution.  This  distribution  is  used  for  the  extension  to  the  3-D  case. 

The  vortex  distribution  for  the  3-D  lifting  surface  as  shown  in  Fig.  2  consists 
of  bound  vortices  and  trailing  vortices.  A  bound  vortex  segment  at  the  ith  chordwise 
strip  and  jth  spanwise  strip  is  a  portion  of  a  U-shaped  vortex  filament,  of  strength 
rtj,  which  propagates  downstream  to  Infinity.  The  wake  tangency  condition  is 
satisfied  at  the  Kutta  point  Kj.  The  vortex  strength  of  the  segment  ij  is 

r 1 j(*i  Sj)  =  T(xlt  sj)  Y (s j)  h(x^ )  1=1,  ...,  1-1 

j  -  1 . J-l  (14) 

where  T  is  the  blade  thickness  distribution,  y  is  the  bound  vortex  strength  in  the 
spanwise  (sj)  direction,  and  h  is  defined  as  the  piecewise  length  of  the  vortex  fil¬ 
ament  along  the  chord.  The  J-l  algebraic  equations  are  formulated  for  y(sj). 

For  lifting  surfaces  without  rotation,  the  trailing  vortex  filaments  are  paral¬ 
lel  to  the  incoming  flow.  For  surfaces  with  rotation,  the  trailing  vortex  filaments 
are  defined  by 

x  =  rt  sinO  ,  y  =  rt  cos0  ,  z  =  zt  +  (15) 

where  the  subscript  t  refers  to  the  trailing  edge,  and  0  is  the  azimuthal  angle 
measured  from  the  trailing  edge. 
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NUMERICAL  ITERATION 


The  two  integral  Eqs.  (10a)  or  (10b)  for  the  source  strengths  and  the  algebraic 
equations  for  the  vortex  strengths  are  coupled  and  should  be  solved  simultaneously. 

In  the  present  approach,  a  global  iterative  procedure  is  used  to  obtain  solutions 
for  both  singularity  strengths.  In  order  to  eliminate  a  large  matrix  inversion,  the 
integral  Eq.  (10)  for  M  is  solved  iteratively  also.  Only  the  small  set  of  algebraic 
equations  for  T  is  solved  by  Gaussian  elimination  procedure.  The  Neumann  iteration 
formula  is  used  in  this  approach  to  solve  Fredholm  integral  equations  of  the  second 
kind.  Hence,  Eqs.  (10a)  and  (10b)  can  be  cast  in  the  form 

4irM(n+1)(P)  -  f(n)(P)  +  Js  [M(n)(Q)K(P,Q)  -  M(n)(P)K(Q,P)]dSQ  (16a) 

for  the  unbounded  flow  case  and 

2irM(n+1)(P)  -  f(n)(P)  +  /s[M(n)(Q)K(P,Q)  -  M(n)(P)K(Q,P)]dSQ 

+  [ 2ir  -  o0(P)]M{n)(P)  (16b) 

for  the  bounded  flow  case,  where  n  is  the  number  of  the  iterations  and  f^n\p)  = 

VN(P)  -  VpN^n\p).  Although  the  iteration  between  sources  and  vortices  is  decoupled 
in  the  global  iteration,  the  interaction  between  them  is  Included  through  the  vari¬ 
ation  of  f(p)  between  iterations. 

One  must  select  a  first  approximation  in  order  to  start  the  iterations.  In 
applying  the  Neumann  iteration  formula,  the  usual  choice  is  M^)(P)  =  VN(P)/2n  for 
Eq.  (6).  However,  as  was  shown  in  [12]  for  the  unbounded-flow  case,  this  choice 
may  result  in  slow  convergence  of  the  iterations  and  a  preferable  first  approximation 
is  given  by  [  1 2 ]  as  H^^(P)  -  V[}(P)/4it.  The  same  value  was  also  selected  for  the 
bounded-flow  case. 

NUMERICAL  EXAMPLES 


Two  sets  of  numerical  results  are  presented  to  illustrate  the  foregoing  proce- 


dure.  The  first  is  for  an  unbounded  flow  case;  namely  a  marine  propeller  tested  In 
open  water  [15,16].  The  second  is  for  a  bounded  flow  case;  namely  an  axial-flow  pump. 

GEOMETRY  INPUTS 


The  surface  S  depicted  in  Fig.  1 ,  on  which  the  sources  are  to  be  distributed,  is 
divided  into  small  panels,  each  with  constant  source  strength.  On  each  panel  a 
representative  point  is  chosen,  where  the  boundary  condition,  Eq.  (1),  will  be 
enforced.  The  panel  area  acts  as  a  weighting  factor  in  the  calculation,  and  the  sur¬ 
face  normal  at  the  representative  point  defines  the  panel  orientation.  The  axisym- 
metric  boundary  surface  is  divided  into  sectors  corresponding  to  the  number  of  blades 
of  the  rotor.  Surface-panel  generation  is  conducted  in  detail  for  only  one  of  these 
sectors.  Panels  for  the  remaining  sectors  are  then  generated  simply  by  rotating  the 
basic  set.  Figures  6  and  7  show  a  total  of  1,650  panels  for  a  3-blade  marine  pro¬ 
peller  in  open  water  and  a  total  of  3,870  panels  for  a  5-blade  axial-flow  pump. 

The  marine  propeller  is  DTNSRDC  controllable-pitch  Propeller  4718,  with  the  diameter 


Fig.  6.  Paneling  on  the  marine  propeller. 
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D  ■  0.61  m,  an  Expanded-Area-Ratlo  (EAR)  of  0.44  and  a  tip  skew  of  20°. 

The  upstream  supporting  strut  used  In  the  measurements  Is  terminated  about  one 
propeller  radius  upstream  of  the  propeller  plane  In  the  calculation  model.  The 
axial-flow  pump  has  a  casing  diameter  D  -  1.04  m.  The  hub  diameter  varies  from 
0.35D  to  0.66D  In  the  blade  section. 


Fig.  7.  Paneling  on  the  axial-flow  pump. 


NUMERICAL  RESULTS 

The  steady  pressure  distribution  was  calculated  on  the  surface  of  Propeller  4718 
at  the  design  advance  coefficient,  .1  =  2nV00/a)D  «  0.751.  Figure  8  shows  the  generated 


Fig.  8.  Wake  vortex  filaments  and  grids  of  one  sector  for 
Propeller  4718  at  J  ■  0.751. 
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wake  vortex  pattern  under  the  operating  condition  by  using  Eq.  (15).  The  wake 
terminates  at  about  1.5  propeller  radii  from  the  trailing  edge.  The  computational 
results  show  that  the  flow  around  the  blade  is  not  affected  by  even  doubling  the  wake 
vortex  filament  length.  Convergence  of  the  solution  was  obtained  after  15  iterations. 
In  Fig.  9,  the  predicted  pressure  coefficients  normalized  with  respect  to  local 
dynamic  pressure  as  defined  in  Eq.  (5)  are  compared  with  experimental  measurements 
[ 15,16]  on  the  blade  surface  at  four  different  radii,  i.e.  2r/D  -  0.5,  0.7,  0.8  and 
0.9,  with  the  propeller  operating  in  uniform  flow.  The  measurements  at  2r/D  =  0.8 
were  only  taken  on  the  suction  side.  Three  pressure  distributions  are  shown  at  each 
radius  of  Fig.  9:  (1)  measured  pressure  distribution,  (2)  distribution  calculated 
by  Kim  and  Kobayashi  [17]  using  a  vortex- lattice  method  (labeled  as  PSF2),  and  (3) 
distribution  calculated  by  the  present  method  (labeled  as  IPF3D).  The  calculated 
results  using  the  present  approach,  i.e.,  IPF3D,  agree  very  well  with  the  measured 
data  except  near  the  trailing  edge  for  2r/D  =0.5  and  0.7.  This  discrepancy  near 
the  trailing  edge  is  due  to  the  wake  tangency  condition  used  at  the  trailing-edge 
Kutta  point,  which  depends  strongly  on  the  direction  given  in  Input  data.  Predic¬ 
tions  of  PSF2  are  lower  than  the  measured  as  well  as  the  current  calculated 
pressures.  In  their  paper,  Kim  and  Kobayashi  also  presented  calculated  and  measured 
data  for  the  propeller  at  two  different  values  of  the  advance  ratio  J.  The  calcu¬ 
lated  and  experimental  results  compare  more  favorably  in  the  case  of  a  NSMB  propel¬ 
ler  than  in  the  case  of  Propeller  4718.  This  is  thought  to  be  due  to  the  following 
reasons:  (l)  the  NSMB  propeller  is  more  heavily  loaded  than  Propeller  4718.  (The 

vortex  lattice  method  [5,6]  usually  predicts  more  accurately  for  a  highly  cambered 
lifting  surface)  and  (il)  Propeller  4718  has  a  tip  skew  angle  of  20°  as  compared 
to  0°  for  the  NSMB  propeller.  (PSF2  uses  induced  velocities  from  other  thickness- 
dominated  components  at  the  propeller  plane  as  inputs  to  represent  the  complete 
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Fig.  9.  Pressure  distributions  on  Propeller  4718  at  J  -  0.751. 


blade  surface  boundary  conditions.  The  deviation  becomes  large  when  the  skew  angle 
Is  large.)  These  two  difficulties  apparently  do  not  occur  in  the  method  presented 
in  this  paper. 


Figure  10  shows  the  suction-side  pressure  distributions  at  two  radii,  l.e. 

2r/D  -0.5  and  0.8,  for  three  off-deslgn  conditions.  Because  the  measurements  were 
only  performed  on  the  suction  side,  the  calculated  distribution  on  the  pressure  side 
Is  not  Included.  Based  on  the  computed  results  on  the  NSMB  propeller  and  a  2-D 
conformal-mapping  calculation,  Hess  [7]  concludes  that  the  leading— edge  pressure 


lb) 

Fig.  10.  Pressure  distributions  on  Propeller  4718 
at  different  values  of  J. 
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Fig.  11.  Wake  vortex  f laments  and  grids  of  one  sector  for 
axial-flow  pump  at  c  -  0.292. 

peak  exists  on  the  pressure  side  for  off-design  condition.  However,  results  from 
the  present  calculation  for  Propeller  4718  do  not  show  this  peak,  which  may  be  due 
to  the  stronger  three-dimensional  relief  for  the  higher  skew  blade  of  Propeller  4718. 

The  axial-flow  pump  used  as  a  calculation  example  is  a  lightly  loaded  machine 
equipped  with  sharp  leading-  and  trailing-edge  blades.  Measurement  data  are  not 
available  for  comparison.  The  calculated  results  are  shown  here  to  demonstrate  the 
feasibility  of  the  methodology  for  an  internal  flow  with  an  axially  rotating  Impeller. 
The  flow  coefficient  c,  defined  as  the  ratio  of  to  the  blade  tip  speed  is  0.292. 
Figure  11  shows  the  wake  pattern  of  the  vortex  filaments,  generated  by  Eq.  (14), 
associated  with  c  ■  0.292.  The  vortex  filaments  extend  1.830  downstream  from  the 
blade  trailing  edge.  The  boundary  condition  for  the  inlet  panels  is  4m/sec.  Figure  12 
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shows  the  computed  pressure  distributions  on  both  suction  and  pressure  sides  of  the 
blade  for  two  strips  parallel  to  the  hub.  The  comparisons  are  made  for  two  differ¬ 
ent  grids  used  in  the  calculations.  The  fine  grid  has  (1580  panels  in  each  blade 
sector)  twice  the  number  of  panels  as  compared  to  the  coarse  grid  (774  panels).  The 
predicted  pressure  distributions,  shown  in  Fig.  12,  for  both  grids  are  obtained 
after  19  global  Iterations.  The  agreement  of  these  results  indicates  the  adequacy 
of  the  coarse  grid  calculation.  Figure  12  also  shows  the  panel  distribution  on  the 
blade  for  the  coarse  grid.  Figure  13  shows  transverse  velocity  distributions,  i.e. 
the  vector  sum  of  the  tangential  and  the  radial  components  of  the  computed  velocity 


(a) 


FRACTION  OF  CHORD 


(b) 


Fig.  12.  Calculated  blade  pressure  distributions  for  axial-flow  pump 
at  c  ■  0.292  for  (a)  5th  strip,  (b)  8th  strip. 
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at  two  sections  of  0.1 ID  and  0.6D  downstream  of  the  trailing  edge.  Due  to  the  pre¬ 
scribed  wake  vortex  filaments  and  the  proximity  of  the  velocity  computing  point  to 
the  vortex  centerline  and  hub  or  shroud  surfaces,  a  few  unrealistically  large  veloc¬ 
ity  vectors  are  seen  In  Fig.  13.  The  mass  flowrates  across  these  sections  are  within 
3-6%  of  the  inlet  condition.  The  calculated  energy  transfer  rate,  defined  as  the 
difference  In  angular  momentum  between  an  Inlet  and  an  outlet,  is  within  8%  of  the 
design  value. 


Fig.  13.  Computed  transverse-velocity  distribution  at  (a)  i  -  0.11D, 
(b)  z  ■  0.60D  from  the  trailing  edge  for  axial-flow  pump. 


CONCLUSIONS 

1.  An  integral-equatlon/slngularity-method  approach  to  the  problem  of  computing  flow 
past  a  rotating  element  with  lift-generating  surfaces  is  feasible  when  using 
source  distributions  for  both  solid  boundaries  and  a  permeable  inlet  flow  control 
boundary,  and  bound  vortex  and  wake  trailing  vortex  filaments  for  lifting 
surfaces. 
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2.  With  a  slight  variation  in  the  governing  integral  equation,  the  method  can  be 
applied  to  both  Internal  and  external  rotating  axial  flows. 

3.  An  effective  mean-line  vortex  strength  distribution,  described  by  Eq.  (14),  has 
been  obtained  for  the  representation  of  the  lifting  surfaces. 

4.  The  Neumann  iteration  scheme,  as  compared  to  the  matrix  Inversion  scheme,  has 
been  shown  to  be  efficient  and  accurate. 

5.  The  global  iteration  between  the  source  strengths  in  the  integral  equations  and 
the  vortex  strengths  in  the  algebraic  equations  converges  fairly  rapidly.  This 
iterative  procedure  was  used  in  conjunction  with  the  Neumann  iteration  scheme  and 
found  to  accelerate  the  convergence  for  both  iteration  processes. 
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Moran 

Feldman 

Lugt 
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Huang 

Nickerson 
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Strieker 

Purnell 

Bein 

Wang 

Henry 

Calvert 
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